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Abstract 

Our proposed method of random phase-free holography using virtual con¬ 
vergence light can obtain large reconstructed images exceeding the size of 
the hologram, without the assistance of random phase. The reconstructed 
images have low-speckle noise in the amplitude and phase-only holograms 
(kinoforms); however, in low-resolution holograms, we obtain a degraded im¬ 
age quality compared to the original image. We propose an iterative random 
phase-free method with virtual convergence light to address this problem. 
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1. Introduction 

Digital holographic display is a promising technique, because a wavefront 
of light scattered from an object can be appropriately reconstructed from the 
display; therefore, this property will enable the realization of an ideal three- 
dimensional display and projector. Holographic projections BBBBB have 
unique properties, including multi-projection B ( ^y which a multi-image is 
projected on multiple screens), projection on screens of arbitrary surface, 
and lensless zoom-able holographic projection (3, S) B- The lensless zoom- 
able holographic projection will lead to the development of an ultra-small 
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projector. Reconstructed images exceeding the hologram size, in general, 
require the random phase; however, this causes considerable problems of 
speckle noise. 

There are well-known methods for improving the random phase applied- 
holograms: the Gerchberg-Saxton (GS) algorithm (ifil. [llj|. the multi-random 
phase method Q. the one-step-phase-retrieval method (OSPR) E3l, and 


pixel separation methods mis ], [i^. The multi-random phase and pixel 
separation methods require display devices with high-speed refresh rates. 
Conversely, random phase-free methods have also been proposed, for exam¬ 


ple, the error diffusion [iJ, ll7|, ll8|, ll9|, l20| and down-sampling methods [21 


These methods can reconstruct clear images; however, the size of the recon¬ 
structed image cannot exceed the size of the hologram, because light fron the 
object does not spread widely j^. Therefore, these methods cannot be used 
for lensless zoom-able holographic projections. 

Recently, new random phase-free methods using virtual special conver¬ 
gence light for amplitude and phase-only holograms (kinoforms) have been 


proposed [22l, 123|, |2J]. Without the assistance of the random phase, this 
method can reconstruct images that exceed the hologram’s size with low- 
speckle noise; however, in low-resolution holograms, we obtain a degraded 
image quality by ringing artifacts, as will be shown in the next section. 

In this paper, we propose an iterative random phase-free method with vir¬ 
tual convergence light to address this problem. Sections 2 and 3 describe the 
proposed method and the results of simulation performed using it. Section 
4 concludes this work. 


2. Proposed method 


m 


The random phase-free method using virtual convergence light is outlined 


22l . l23l . l24| . The calculation setup is shown in Fig. [H This method is 


applicable to amplitude computer-generated holograms (CGHs) and kino¬ 
forms. Instead of using the random phase, the complex amplitude on the 
image plane Ui(xi,yi) is multiplied using virtual convergence light given by 


w{xi,yi) = exp(-i7r(a:- -F y1)/\fi), 


( 1 ) 


where fi = zi + is the focal length, z\ is the distance between the focus 
point of the convergence light and the hologram, and is set to the distance 
at which the hologram just fits to the cone of the convergence light, and Z 2 
is the distance between the object and hologram. 
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Figure 1: Calculation setup for the random phase-free method using virtual convergence 
light. 


Here we describe how to determine /^. Using a simple geometric relation, 
we can derive Sh/'2 \ Si/2 = Zi \ fi where the areas of the image and the 
amplitude CGH (or kinoform) are given by Si x Si and ShX Sh^ respectively. 
Therefore, we obtain: 

fi^Z2/il-Sh/Si). (2) 

To avoid overlap between the reconstructed image and the 0-th order light, 
the original object must be shifted from the optical axis by a distance of o. 
Owing to the addition of this shift amount, the focal length of the convergence 
light is expressed as 

fi = Z 2 I (1 — Shj (5'i + 2o)). (3) 

We calculate the complex amplitude in the hologram plane by 

Uh{xh, Vh) = Pvop^^{ui{xi, yi)w{xi, y^)}, (4) 

where Propz 2 {-} denotes the diffraction calculated at the propagation dis¬ 
tance, Z 2 . The following equation is used to calculate the amplitude CGH, 

I{xh,yh) from Uh{xh,yh)- 


I{xh,yh) = ^{uhixh,yh)}, 


( 5 ) 
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where 3?{-} denotes the real part of Uh{xh,yh). In calculating the kinoform, 
6(xh,yh), from Uh{xh,yh), the following equation is used: 

0{xh, Vh) = fMg{uh{xh, yh)}, (6) 


where arg{-} indicate the use of only the argument of the complex ampli¬ 
tude. Unfortunately, in kinoforms, the reconstructed images using the ran¬ 
dom phase-free method are heavily degraded. This is a common problem of 
kinoforms 0,[2l|. To overcome this problem, we apply the error diffusion 
method dal Enl | to the random phase-free kinoform as follows: 


0{xh, yh) = ^'Q{0{xh, yh)}, (7) 


where ED{-} denotes the error diffusion operator. The details of the error 
diffusion are described in Section 12.11 For the details of the random phase 
free-amplitude CGH and kinoform, refer Ref. 


22, 23 


Figure (21 shows reconstructed images from the amplitude CGH and kino¬ 
form. Figures[2](a) and[2](c) are the reconstructions from the amplitude CGHs 
of the images “Lena” and “Mandrill’, respectively. Figures [21(b) and [21(d) are 
the reconstructions from the kinoforms of the images, respectively. These 
reconstructed images do not have the speckle noise induced by the random 
phase; however, ringing artifacts overlap in these reconstructed images. 


2.1. Improved image quality using an iterative method 

To overcome the problem with ringing artifacts problem in the recon¬ 
structed images, we apply an iterative enhancement of the random phase-free 
method. Figures El and 01 show the improved method for amplitude CGHs and 
kinoforms using the iterative algorithm, respectively. The kinoform yielded 
by the iteration method is almost the same as that in Fig. 01 except for the 
error diffusion. The calculation steps are as follows: 

1. We start the iteration by multiplying the virtual convergence light 
w{xi,yi) by the original image ai(xi,yi). The image plane is denoted 
by Ui{xi,yi) = ai{xi,yi)w{xi,yi). 

2. We calculate the scaled and shifted diffraction, which can be used to 
obtain light propagation at different sampling pitches on the image 
and hologram, as well as the off-axis propagation, by Uh{xh,yh) = 
Prop^.,I Uj( Xj.U i)wiXj.Vi)}. The diffraction we used is ARSS Fresnel 
diffraction (25[, which is a scaled and shifted diffraction. The sampling 
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Figure 2: Reconstructed images from amplitude CGH and kinoform. (a) and (c) recon¬ 
structed images from amplitude CGHs: (b) and (d) reconstructed images from kinoform: 
These reconstructed images are contaminated by ringing artifacts. 


pitches on the image and hologram planes are set to pi and ph^ respec¬ 
tively. The amount of shift away from the optical axis on the image 
plane is {o^,Oy).] 

3. We extract only the real part of Uh(xh,yh) for the amplitude CGH 
or only the argument of Uh(xh,yh) for the kinoform (which represents 
constraints in the hologram plane). 


9{xh,yh) 


^{uh{xh,yh)} (For amplitude CGH) (8a) 

ED{avg{uh{xh,yh)}} (For kinoform) (8b) 


4. We back-propagate Uh{xh, yh) to the image plane, u[{xi, yi) = VYo\)_^^{g{xh, yh)} 
by ARSS Fresnel diffraction. If we use the offset parameter, the shift 
amount away from the optical axis should be set to (— Oa,, —Oy). 
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Figure 3: Iterative method for amplitude CGH. 

5. We update only the amplitude of u[{xi^yi) according to 

Ui{xi,yi) = ai{xi,yi)u[{xi,yi)l\u[{xi,yi)\. (9) 

This is the constraint in the image plane. Subsequently, we obtain the 
new image plane Ui{xi,yi). 

We repeat the above iteration steps (2) to (5) until reaching the preset iter¬ 
ation number or a good image quality. 

In the step 3, the Floyd-Steinberg error diffusion method we used was 
calculated using 


0{xh,yh + l) ^ 6{xh,yh + 'i-) + wie{xh,yh), ( 10 ) 

9(xh + l,yh — l) d{xh + l,yh — 1) + W2e{xh,yh), ( 11 ) 

9{xh + l,yh) ^ 9{xh + l,yh) + W3e{xh,yh), ( 12 ) 

9{xh + l,yh + l) 9{xh + l,yh + 1) + W4e(xh,yh), (13) 


where the coefficients are rci=7/16, rc2=3/16, W3=b/1Q, and ?n4=l/16 (l^ : 

indicates that the current value is overwritten, and e{xh, yh) = Uh{xh, yn) — 
9{xh,yh)- Before applying the error diffusion, we may need to normalize 
the complex amplitude Uh{xh,yh) by Uhixh.yh) ^ Uh{xh,yh)/d, where d = 
max.\uh{xh,yh)\ is the maximum absolute value in the kinoform. 
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Figure 4: Iterative method for kinoform. 


3. Results 

We show the effectiveness of the proposed method through simulation. 
The calculation is performed under conditions where the resolution of the 
images is 2, 048 x 2, 048 pixels, the wavelength is 532 nm, the pixel pitches 
on the image plane and the hologram plane are p* = 8/um and ph = 4pm, 
respectively, the distance between the image and hologram planes is 1 m, and 
the amount of shift away from the optical axis is Oy) — (20mm, 20mm). 

The upper images in Fig. O show reconstructions of the image “Lena” 
from amplitude CGHs, and the lower images show reconstructions kinoforms, 
which are generated by the proposed method. Increasing the iteration num¬ 
ber results in improvement of the image quality in both amplitude CGH and 
kinoform. Figure IHl shows the peak signal-to-noise ratio (PNSR) between the 
original image and the reconstructed image when increasing the number of 
iteration. The blue and green lines indicate the PSNRs for the amplitude 
CGHs and kinoforms, respectively. The improvement effect converges at five 
iterations, and the PSNRs of the amplitude CGH and kinoform are about 30 
dB and 27 dB, respectively. 

Figure [7| shows the reconstructions of the image “Mandrill”. The upper 
and lower images are obtained from amplitude CGHs and kinoforms gen¬ 
erated by the proposed method, respectively. Figure [HI shows the PNSR 
between the original and the reconstructed images when increasing the num¬ 
ber of iteration. The improvement effect is the same as that in FigjHl and the 
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PSNRs at five iterations of the amplitude CGH and kinoform are about 23 
dB and 22 dB, respectively. Overall, the proposed method greatly improves 
the image quality. 




Reconstructed images from amplitude CGHs 


Reconstructed images from amplitude kinoforms 


Figure 5; Reconstructions of the image “Lena”. The upper and lower images are recon¬ 
structed from amplitude CGHs and kinoforms, respectively. They are generated by the 
proposed method. 

Figure ini shows the reconstructed images from amplitude CGHs with ran¬ 
dom phase (left) and with the proposed method (right) at five iterations. The 
PSNRs for the left and right images are 9.8 dB and 30.5 dB, respectively. 
Figure [ini shows the reconstructed images from kinoforms with random phase 
(left), with the proposed method not applying the error diffusion (middle) 
and with the proposed method using the error diffusion (right) at five itera¬ 
tions. The PSNRs for the left and right images are 9.0 dB, 9.4 dB, and 27.1 
dB, respectively. 

4. Conclusion 

We proposed an iterative technique for generating amplitude CGHs and 
kinoforms using the random phase-free method with virtual convergence 
light. The proposed method succeeded in mitigating ringing artifacts in the 
reconstructed images. 
























Figure 6: PSNRs for the reconstructed images of “Lena” from the amplitude CGHs and 
kinoforms. 
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